Slow imbalance relaxation and thermoelectric transport in graphene 
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We compute the electronic component (k) of the thermal conductivity and the thermoelectric 
power (a) of monolayer graphene, within the hydrodynamic regime, taking into account the slow rate 
of carrier population imbalance relaxation. Interband electron-hole generation and recombination 
processes are inefficient due to the non-decaying nature of the relativistic energy spectrum. As a 
result, a population imbalance of the conduction and valence bands is generically induced upon the 
application of a thermal gradient. We show that the thermoelectric response of a graphene monolayer 
depends upon the ratio of the sample length to an intrinsic length scale Iq, set by the imbalance 
relaxation rate. At the same time, we incorporate the crucial influence of the metallic contacts 
required for the thermopower measurement (under open circuit boundary conditions), since carrier 
exchange with the contacts also rela^xes the imbalance. These effects are especially pronounced for 
clean graphene, where the thermoelectric transport is limited exclusively by intercarrier collisions. 
For specimens shorter than Iq, the population imbalance extends throughout the sample; k, and a 
asymptote toward their zero imbalance relaxation limits. In the opposite limit of a graphene slab 
longer than Iq, at non-zero doping k and a approach intrinsic values characteristic of the infinite 
imbalance relaxation limit. Samples of intermediate (long) length in the doped (undoped) case are 
predicted to exhibit an inhomogeneous temperature profile, whilst k and a grow linearly with the 
system size. In all cases except for the shortest devices, we develop a picture of bulk electron and hole 
number currents that flow between thermally conductive leads, where steady-state recombination 
and generation processes relax the accumulating imbalance. Our analysis incorporates, in addition, 
the effects of (weak) quenched disorder. 

PACS numbers: 



I. INTRODUCTION 

Both quenched disorder and interparticle interac- 
tion effects influence electric and thermal transport in 
graphene J^^^^^di^n^nL^ i ^ o i ^ ^ At exactly zero doping (the 
so-called "Dirac point" ) , the electron-hole fluid of mass- 
less Dirac quasiparticles is predicted to exhibit a perfectly 
finite, non-zero dc electrical conductivity a for temper- 
atures T > 0, even in the absence of disorder, due en- 
tirely to electron-hole collisions i^'^i^°'^^i^^ Since the car- 
rier plasma is electrically neutral at zero doping, an ap- 
plied electric field does not couple to the centcr-of-mass 
momentum of the fiuid; instead, electrons and holes are 
driven in opposite directions, and electron-hole collisions 
limit the developing electric current. By contrast, a tem- 
perature gradient can induce a thermal drift of both elec- 
tron and hole fluid components in the same direction; 
in the absence of quenched disorder, the resulting en- 
ergy current grows unimpeded by interparticle collisions, 
which cannot influence the center of mass momentum. 
It has been therefore claimed that the electronic com- 
ponent K of the thermal conductivity of clean, undoped 
graphene is infinite^ being limited only by the amount 
disorder in a dirty sample. We note that the thermoelec- 
tric power a = at the Dirac point due to (approximate) 
particle-hole symmetry. 

In this paper we demonstrate that the above picture is 
incomplete: a clean graphene sample at the Dirac point 
will in general exhibit a perfectly finite electronic ther- 
mal conductance Gtu, independent of the sample length 
L for L ^ Iq, where Iq represents a certain intrinsic 



length scale (set by the bulk properties of the material 
and the average temperature). In such large specimens, 
the response is predicted to be inhomogeneous, with tem- 
perature gradients confined to boundary regions of size 
Iq adjoining the sample edges. In the opposite limit 
L Iq, the temperature falls linearly across the en- 
tire sample, and we find a well-defined, L-independent k. 
At non-zero doping, the thermoelectric transport prop- 
erties exhibit several different behavioral regimes. In the 
large system size limit L — > oo, both n and a asymp- 
tote toward perfectly finite values at non-zero T and ^, 
even in the absence of disorder, in accord with the re- 
sults of Ref. [l^. In the opposite limit of a sufficiently 
short device L <C Lq, we obtain new, completely differ- 
ent results. The mechanism responsible involves only in- 
terparticle collisions, but in particular the comparatively 
slow relaxation of electron-hole population imbalance. 

Imbalance refers to a state in which the electron and 
hole populations deviate from their values in chemical 
equilibrium, so that these carriers possess independent 
chemical potentials ^e,h and particle densities n^^h- By 
contrast, equilibrium slaves /.te = — /x/j = y^t in graphene 
(a zero bandgap semiconductor), so that both rig and Uh 
are completely determined by /i and the temperature T. 
Recombination or generation processes relax a popula- 
tion imbalance within a time interval tq, the imbalance 
relaxation lifetime. The lowest order ("Auger" or two 
particle collision) relaxation processes typically involve 
the absorption or emission of an electron-hole pair, i.e. 
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FIG. 1: Kinematical constraints for imbalance relaxation. 
The left figure shows the Feynman diagram for the typi- 
cal two-particle decay process given by Eq. (|l.lap : Si and 
Pi respectively denote the energy and momentum of the i**^ 
electron or hole. The right figure depicts momentum con- 
servation for this process. The length of the dashed path 
Lf = IpjI + \pg \ + IP4I traced out by the sum of "decay prod- 
uct" momenta is always greater than or equal to the length 
Li = IpjI of the "parent" particle momentum: Lf > Li. If the 
particle energy spectrum takes the form e(p) = |p|'', then the 
depicted decay process is kinematically forbidden for /3 < 1. 
For /3 = 1, only forward-scattering is allowed. 
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In graphcnc, however, these processes are kinematically 
suppressed by the conservation of energy e and momen- 
tum p, because the spectrum e{p) of the quasiparticles is 
not decaying 
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{Negative curvature of the spectrum at T = 
K arises due to the logarithmic renormalization of 
the Fermi velocity vp attributed to electron-electron 
As explicated in Fig. [U the linear 
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spectrum allows only decay products with coUinear mo- 
menta (i.e. pure forward scattering), but these processes 
make a negligible contribution to the imbalance relax- 
ation. The sublinear spectrum forbids even this forward 
scattering decay. In clean graphene, higher order (e.g. 
three particle collision) imbalance relaxation processes 
are already allowed, while impurity-assisted collisions will 
contribute in a disordered sample; tq is therefore likely 
finite, although it may significantly exceed other relax- 
ation times in the system. 

In the limit of zero relaxation, a graphene monolayer 
probed through thermally conducting, electrically insu- 
lating contacts would possess electron and hole popu- 
lations that are strictly conserved. In direct analogy 
with a single component, non-relativistic classical gas^^ 
the electron-hole plasma in clean graphene with vanish- 
ing imbalance relaxation would exhibit a finite electronic 
thermal conductivity k for arbitrarily large system sizes. 
(The imbalance relaxation length Iq, introduced above, 
diverges as tq — > 00.) In this regime, interparticle colli- 
sions facilitate heat conduction without particle number 
convection. In this paper, we will demonstrate that the 
same behavior obtains for non-zero imbalance relaxation 
{1/Iq > 0) in the limit of short samples, L <^ Iq. By 
comparison, prior workiS effectively assumed infinite re- 



laxation of population imbalance {Iq — > 0). We demon- 
strate that the results previously obtained in Ref. [l^ for 
both K and the thermoelectric power a at non-zero dop- 
ing and temperature emerge in the limit of asymptoti- 
cally large system sizes L ^ Iq for a finite rate of im- 
balance relaxation, {Iq > 0). We will show that when 
the ends of a graphene slab of length L > Iq are held at 
disparate temperatures and no electric current is permit- 
ted to flow, steady state particle convection does nev- 
ertheless occur; carrier flux is created or destroyed by 
imbalance relaxation processes near the terminals of the 
device. Thermopower measurements require the junc- 
tion of the graphene slab with metallic contacts. We 
incorporate into our calculations carrier exchange with 
non-ideal contacts, which also relax the imbalance, and 
we carefully delimit the regimes in which deviations from 
the infinite relaxation limit should be observable in ex- 
periments. The effects of weak quenched disorder are 
included in all of our computations. 

In this paper, we restrict our attention to the hy- 
drodynamic (or "interaction-limited" ) transport regime, 
where inelastic interparticle collisions dominate over elas- 
tic impurity scattering. Prior work addressing ther- 
moelectric transport in the opposite, "disorder-limited" 
regime, in which real carrier-carrier scatte ring pr ocesses 
may be neglected, includes that of Refs. IlllllTlflsl . In 
the disorder-limited and a are determined by 

the energy-dependence of the electrical conductivity, via 
the "generalized" Wiedemann-Franz law and Mott re- 
lation, respectivelyfi^ The effect of the slow imbalance 
relaxation upon the dc conductivity in graphene under 
non-equilibrium interband photoexcitation has also been 
addressed/^ 

What essential new physics emerges through the in- 
corporation of imbalance relaxation effects into the de- 
scription of thermoelectric transport, and how can it 
be extracted from experiments? The entirety of linear 
transport phenomena in graphene within the hydrody- 
namic regime is essentially quantified by four intrinsic 
parameters. A finite rate of imbalance relaxation means 
that electrons and holes respond independently to ex- 
ternal forces; the single "quantum critical" conductiv- 
ity identified previously in Refs. [sHolfTol generalizes to a 
2x2 tensor of coefficients, with diagonal elements CTee, 
ahh and off-diagonal elements a eh = cThe] all are me- 
diated entirely by inelastic interparticle collisions. The 
description is similar to that of Coulomb dragi^ the di- 
agonal element dee {<^hh) characterizes the response of 
the conduction band electrons (valence band holes) to a 
(gedanken) electric field that couples only to that carrier 
type, whereas a^h characterizes the "drag" exerted by 
one carrier species upon the other (due to electron-hole 
collisions) under the application of such a field, a^e and 
ahh are related by particle-hole symmetry, leaving two 
independent parameters which we can take as Ce/i and 
(Tmin = (cee +o'hh ^ '^<^eh)] the latter Combination is equal 
to the bulk dc electrical conductivity a at the Dirac point 
in the hydrodynamic regimcji^ The imbalance relaxation 
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TABLE I: Intrinsic transport parameters for graphene in the hydrodynamic regime, and how to extract them from ther- 
moelectric transport (and other) measurements. Hi{L) is the electronic contribution to the thermal conductivity, and a{L) is 
the thermopower. Here and below, L denotes the length of the putative graphene device, which should be compared to the 
imbalance relaxation length Iq. The parameters Iq and a^h have been introduced in this work. Relevant results obtained in 
this paper can be found in the equations listed in the third column. 



Parameter Description 

Iq Imbalance relaxation length 

Zei Elastic mean free path 

(due to quenched disorder) 

(Tmin Inelastic carrier-carrier scattering 

= (Tee + (Thh — 2aeh coutribution to the bulk a 

(Teh Off-diagonal ("drag") component 

of the conductivity tensor due to 
inelastic electron-hole scattering 



Extract from measured quantity 

• Temperature T{x) or electrochemical potential V'S>{x) profile 
[Eqs. (1133, (HSU, (IOD1) . and (HHJ] 

• Bulk dc conductivity a [Eq. p. 211) ] , or 

• /t(i) — > i^aa and oi{L) — > Qoo obtained in the limit Iq 
[Eqs. fl?^ and 

• Minimum conductivity at the Dirac point [Eq. (|2.21|l ]. or 

• k(I/) — » Kao and a{L) — > aaa obtained in the limit Iq 
[Eqs. ip:^ and 

• ft{L) or a{L) obtained in the limit L <^ Iq 
[Eqs. (|43j and gH)] 



length Iq constitutes the third intrinsic graphene param- 
eter, while the fourth is provided by the elastic mean free 
path Zei for a disordered sample. 

Of these, le\ and (Tmin can both be obtained from either 
the bulk dc conductivity a (measured at variable doping) , 
or from the combined measurement of the "infinite imbal- 
ance relaxation" thermoelectric transport coefficients Kqo 
and aoo- Here, Kqo and ctoo are the (electronic) thermal 
conductivity and thermopower, respectively, as measured 
in a graphene device with L ^ Iq. By contrast, the off- 
diagonal parameter cTeh can be independently determined 
through measurement of either k or a, for a sufficiently 
short device satisfying L <^ Iq. (A precise discussion of 
the limiting and crossover behaviors of k and a, incorpo- 
rating the complicating effects of the external contacts, 
is presented in Sec. IIVI of this paper). Finally, the im- 
balance relaxation length Iq can be ascertained via the 
measurement of either k or a in the crossover regime 
L ^ Iq, or through spatial resolution of the inhomoge- 
neous temperature or electrochemical potential profiles 
across a device with L > Iq. The four intrinsic graphene 
transport parameters discussed above and suggestions for 
their experimental determination are listed in Table U 

The rest of this paper is organized as follows. 

In Sec. ini we formulate a hydrodynamic description 
of carrier transport that admits carrier population im- 
balance and relaxation. In order to illustrate ideas, 
we apply this formalism to a putative experimental de- 
vice. The hydrodynamic approach allows us to obtain 
the inhomogeneous temperature, chemical potential, and 
number current profiles across the device; some of these 
are sketched in Fig. |4] for a system at zero doping. In 
Sec. mil we derive the thermal conductance Gth at the 
Dirac point, and we discuss the short and long device 
asymptotics of heat transport. General results for k and 



a at arbitrary doping are obtained and discussed in the 
penultimate Sec. HVl 

We neglect phonon effects throughout the body of 
this work, but these are considered in the Conclusion 
Sec. |Vl We focus in particular upon the infiuence of 
electron-phonon interactions on the electronic transport 
coefficients derived in this paper. We demonstrate that 
phonons may be neglected at low temperatures or for 
graphene samples of mcsoscopic size. 



II. HYDRODYNAMIC FORMULATION AND 
SOLUTION 

A. Two fluid hydrodynamics 

We take as our starting point the hydrodynamic equa- 
tions of motion for carriers in graphene, expressed in rel- 
ativistically covariant notation: 

d.Je = ~dai - -el, (2.1) 

where J* and denote the electron and hole elec- 
tric 3-current densities, respectively, G'-' is the trace- 
less energy-momentum tensor for the (classically) scale- 
invariant two-component plasma, and F^^ is the Faraday 
tensor, incorporating both external and self-consistent 
fields. The Fermi velocity in graphene is denoted by 
vp, while e = — jej is the electron charge. Summation 
over repeated "spacetime" indices is implied in Eqs. (|2.1[) 
and (|2.2p . where G {x°,x^,x'^} = {vFt,x,y}, while 
di = d/dx\ The quantity I in Eq. (|2.ip is the imbal- 
ance relaxation flow, which is proportional to the rate 
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1 Itq of carrier recombination or generation between the 
electron and hole bands. The frictional force density 6*, 
in Eq. (|2.2p manifests the effects of elastic scattering of 
carriers by quenched disorder. In a microscopic quantum 
kinetic equation treatment, / and obtain as certain 
momentum averages of the inelastic and elastic collision 
integrals, respectively. In the limit / — > 0, electrons and 
holes are separately conserved. 

By adopting the hydrodynamic approach, we posit fast 
equilibration of the electron-hole plasma due to inelas- 
tic collisions; in particular, we assume Tj^ < tv i (the 
"interaction-limited" transport regime, see Ref . KLlh . with 
Tin the inelastic lifetime due to interparticle scattering, 
and Tel the elastic lifetime due to quenched disorder. 
We express J\ and O*-' in terms of local thermody- 
namic variables and a hydrodynamic 3-velocity [/' = 
7(?;f,u), with u the ordinary fluid velocity and 7^ — 
(1 — u^/u|,)"^. Incorporating dissipative deviations from 
local equilibrium, we write^ 

Jl = e (neC/* + i^:) , (2.3) 
Jl = + J^;,) , (2.4) 



9tp + V • J = 0. 



(2.9) 



(2.5) 



where v\ ^ and Q"^^ define dissipative fluctuations of the 
electron and hole number currents and stress tensor, re- 
spectively. In Eqs. (|2.3p and (|2.4p . rie and nn represent 
proper (rest frame) electron and hole densities. V de- 
notes the total pressure in Eq. (|2.5p , where we have used 
the thermodynamic relation 37-" = (), with f) the enthalpy 
density; this is a consequence of relativistic scale invari- 
ance. In this same equation, g*-' denotes the Minkowski 
metric tensor 

In what follows, we will neglect for compactness of pre- 
sentation the non-equilibrium component of the stress 
tensor , which describes viscous effects in non-uniform 
fluid flow. For the thermoelectric transport problem 
studied here, viscous drag along the sample edges is ir- 
relevant under the condition W"^ ^ riVpTs\/3V, where 
W is the sample width transverse to parallel electric and 
heat current flows, and rj is the first (dynamic) viscosity 
coefficient.— (The second viscosity ( vanishes for a mass- 
less relativistic gas.)i^ By inserting the decomposition in 
Eqs. ([01) " ([2?5l) into Eqs. ([2?T|) and ([2?2l) . and taking the 



non-relativistic limit u ^ Vp, we derive the following 
entropy, momentum, and particle number balance equa- 
tions: 



T 



dtS+V ■ ysU - —Ue - —h>h 

Vh VL- bel , 



T 



3Pdu „ „ „ , 

^ — = pE - rig Vpe - rihVfih - s vT + bei, 
Vp at 

dm + V -Jn^ ~2I, 



(2.6) 

(2.7) 
(2.8) 



Here, s is the entropy density, fie.h are the electron and 
hole chemical potentials, and E is the electric field. J„ 
and J respectively denote the total carrier number and 
electric current densities in Eqs. (j2.8p and ()2.9|) : these are 
defined as 



nu + Ve 

Je + J/i 
pu + e{Ve 



In Eqs. (l^- dmil . 



rih 



p = e{ne - Uh) 



(2.10) 
(2.11) 

(2.12) 



represent the net carrier number and electric charge den- 
sities, respectively. Let us also define the imbalance fii 
and relative p, chemical potentials, via 



Pi 



Pe + Ph 



p 



Pe - Ph 



On the left-hand side of Eq. (^7)) . 



- = 9t + u ■ V 

at 



(2.13) 



(2.14) 



is the material derivative. 

Within linear response, Eq. (|2.6|) implies that 
the thermodynamic "forces" {eE - TV(pe/T), -eE - 
T'V{ph/T),u,2pj} determine the conjugate "fluxes" 
{fe^/i, bel, /} via the matrix equation 



"9 






M 


bel 




/ 





T J 

TV 



(2.15) 



eE-TV 
-eE 

u 

2pi 

Onsagcr reciprocity^^ of the entropy balance Eq. 
dictates that the kinetic coefficient matrix M is symmet- 
ric, 

M"^ = M, 

where the superscript "T" denotes the matrix transpose 
operation. Assuming vanishing or weak quenched dis- 
order and slow imbalance relaxation, M can be written 
as 



M 



Oeh CThh 























V'j. Tel 









2h -1 



(2.16) 



The electric conductivities (Tab, ab £ {ee, hh, eh}, arise 
solely due to interparticle collisions, and can be com- 
puted in principle within a microscopic quantum kinetic 
equation (QKE) formulationji^, The elastic lifetime Tei de- 
termines the frictional force density bei oc — u; Eqs. (|2.15p 
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and (j2.16p provide an implicit definition for Tei. Finally, 
Xq is a dimensionless parameter that characterizes the 
efficacy of generation and recombination processes. 

Both {(Tab} and Xq are functions of the dimensionless 
ratios ^./ksT and fii/ksT. Particle-hole symmetry re- 
quires that the diagonal conductivity elements satisfy the 
condition 



ficient matrix in Eq. (|2.16p . 



JJ_ P^i \ 



ksT ksT J 



(2.17) 



Eqs. (|2.15|) and (|2.16p therefore assert that the interpar- 
ticle collisions give rise to two independent kinetic co- 
efficients (Tee and Geh that carry the units of electrical 
conductance, in addition to the dimensionless imbalance 
relaxation parameter Xq. In this paper we focus upon 
transport in the non-degenerate regime, fc^T 3> 
within the accuracy of the linear response approximation, 
{(Tab} and Ag can then be regarded as fixed constants, 
typically evaluated at the Dirac point (fi — ~ 0). Un- 
der these conditions, Eq. (j2.17p implies that a^e = crhh- 

It is worth pointing out that the friction density bgi and 
imbalance relaxation flow / appear already at the level 
of the "ideal" hydrodynamics: to compute (in principle) 
the parameters 1/rei and Xq in Eq. (|2.16p . it is suffi- 
cient to solve the associated QKE at zeroth order in the 
inelastic relaxation time Tin, with electron and hole dis- 
tribution functions locally constrained to take the Fermi- 
Dirac formj ^^i^'^i^^ The ideal hydrodynamics is captured 
by Eqs. ([^ - ([m]) with Ve = i^h 0. 

To zeroth order in t\„, bei is nonvanishing for a con- 
vective particle flow in the rest frame of the disorder 
(u ^ 0), while nonzero / arises whenever a population 
imbalance occurs (/z/ ^ 0). The explanation for this is 
as follows: the inelastic collision integral governing in- 
terparticlc scattering for a clean graphcne system with 
vanishing imbalance relaxation would possess zero modes 
associated to homogeneous fluid convection (momentum 
conservation) and global shifts of the total carrier number 
density. These zero modes are made "massive" through 
the introduction of disorder and the inclusion of some 
mechanism for imbalance relaxation (such as three parti- 
cle collisions). For graphcne in the hydrodynamic regime, 
the "masses" (scattering rates) (l/rei) and (I/tq) respec- 
tively associated to weak disorder and inefficient imbal- 
ance relaxation are much smaller than inelastic scattering 
rate {1/t\„) for electron and hole number-conserving col- 
lisions. Even in the limit of arbitrarily efficient equilibra- 
tion due to frequent and strongly inelastic such collisions 
(Tin — 0), nonvanishing u and /i/ induce dissipation by 
coupling to these weakly massive modes. The dominant 
effect of these forces is the generation of bei and / in the 
ideal hydrodynamic description. 

By contrast, a small but non-zero t\„ allows out-of- 
equilibrium deformations of the electron and hole dis- 
tribution function shapes. The {(Tab} acquire nonzero 
values in the first order of the QKE expansion in t\„. 
This separation of zeroth and first order responses justi- 
fies the assumed block-diagonal form of the kinetic coef- 



B. 



"Bulk" kinetic coefficients 



In the limit of infinite imbalance relaxation, Xq oo 
in Eq. (j2.16p . we must slave fi^ = —fJ-h = M and ^/ = in 
Eqs. (|2.6p . p.7p . and ()2.15p . Assuming steady-state con- 
ditions, Eqs. (l2Jl)~([Ml), jSHSl), and ((2J6l) allow the com- 
putation of the electric J and heat Jg current densities 
in the presence of electrochemical potential and tempera- 
ture gradients. The heat current implied by the left-hand 
side of the entropy balance Eq. (j2.6p is 



(2.18) 



J, = TSU - I^LeVe - f^hl^h 
= SVU - fljJn - -J, 



where on the second line we have used Eqs. (|2.10p and 
p. lip . The linear response takes the usual form2£ 



J = (TE + aa^ (- VT) , 



(2.19) 



^ J, = aa^e + + '^c^'L ) (- VT) , (2.20) 



where 



e = E - -Vyu 



is the electrochemical potential gradient. The thermo- 
electric response in Eqs. (|2.19p and (|2.20p is characterized 
by the bulk dc electrical conductivity a, and the infinite 
imbalance relaxation limits of the thermopower aoo and 
thermal conductivity Kqo ■ In terms of the intrinsic kinetic 
parameters defined via Eq. (|2.16p . one finds that 



(Tmin I 2^ 
3"P«F^e|(Tmm 

Ta 

vpldP _ JJ_ 
Ta eT' 



(2.21) 
(2.22) 
(2.23) 



where the minimum conductivity at the Dirac point is 
given by 



(T/ih — 2a, 



eh ■ 



(2.24) 



In these equations, Zgi = vpT^] is the elastic mean free 
path due to quenched impurity scattering. In the clean 
limit, the thermopower in Eq. (|2.23p simplifies to aoo = 
s /p, which may be interpreted as the "transport entropy" 
per charge.— The results of Eqs. (|2.2ip - p.23p were orig- 
inally obtained in Ref. [l^- 

All three thermoelectric coefficients in Eqs. (|2.2ip - 
([2:23)) depend only upon the combination (a^e + cFhh — 
2(Te/i), for arbitrary doping. We will demonstrate be- 
low that these same formulae apply in the presence of 
finite imbalance relaxation, in the limit of a sufficiently 
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T = T, 



V(0) 



V(L) 



J 



FIG. 2: Schematic setup for the two terminal thermopower 
and thermal conductivity measurements. The graphene slab 
of length L is denoted by the dashed gray line. The contacts 
at a; = 0, L are respectively held at temperatures T — Ti,T2 
via external thermostats. The two contacts are electrically 
bridged by a highly resistive wire, represented in the figure 
by the thick black interconnect. The development of a ther- 
moelectric response V{L) — V{0) 7^ in the graphene induces 
a voltage drop AT4, = — in the wire, which drives 
the current J. Measurement of J (via galvanometer) and 
knowledge of the wire resistance allows computation of the 
graphene thermopower, which is simply the ratio of AVu, to 
the temperature drop AT = Ti — T2, in the limit of infinite 
wire resistance. 



long graphene sample. In the opposite limit of a short 
graphene device, we obtain new, completely different re- 
sults for both a and , clS detailed in Sec. IIVI Combined 
with particle- hole symmetry [Eq. (|2.17p ] , measurement of 
long and short graphene devices should allow an exper- 
imental determination of all three parameters (Tee, <^hh, 
and aeh- 



C. Experimental geometry 

We consider an experiment in which a rectangular strip 
of graphene is terminated with metallic contacts at op- 
posite ends of the strip; we take the strip to lie along 
the X-axis between x = and x = L. In a thermal con- 
ductivity measurement, the leads are held at different 
temperatures, r(0) = Ti and T{L) = T2, and the heat 
current 3g [Eq. (|2.18[) ] is measured. For convenience, let 
us introduce 



T 



T1+T2 



AT = Ti- T2 



(2.25) 



To determine the thermopower, we analyze a gedanken 
measurement in which the metallic contacts are intercon- 
nected by a highly resistive "wire" ; we sketch a schematic 
setup in Fig. [2] The wire is taken to be composed of a 
disordered metal that is approximately particle-hole sym- 
metric, and thus manifests no thermoelectric voltage of 
its own. Let V{x) denote the electric potential profile 
within the graphene. The application of the tempera- 
ture gradient AT 7^ induces a drop V{L) — V{0) across 
the graphene slab via its thermoelectric effect. Through 
local electrochemical quasiequilibration at the contacts. 



the graphene potential drop translates into a voltage dif- 
ference AVyj = — across the ends of the wire. 
Here, V^^ and denote the electric potentials near the 
contacts situated at a; = L and a; = 0, respectively, just 
inside of the metal wire (outside of the graphene) . A 
drives an electric current through the wire linking the 
contacts. We assume local clcctroneutrality throughout 
the wire; therefore, the diffusion component of the elec- 
tric current vanishes outside of the graphene. In the limit 
of arbitrarily large wire resistance, the thermopower is 
then simply the ratio 



a 



AK, 
AT ■ 



(2.26) 



Knowledge of the wire resistance and a galvanic measure- 
ment of the current thus allow experimental determina- 
tion of a. 

We specialize to the quasi- ID strip geometry discussed 
above, and assume a steady-state linear response to the 
application of unequal temperatures at the contacts. A 
small AT induces proportional deviations from zero of u, 
/i/, J„, and V{x); we define the graphene electric poten- 
tial such that V{x) ^0 in the limit AT ^ 0. The tem- 
perature and relative chemical potential are expanded 
about their average values: 

T{x)=T + -^{x), ^i{x) = fi + <i>{x), (2.27) 

where and <I> are assumed proportional to AT. 



D. Boundary conditions 

We impose the following boundary conditionsS^ relat- 
ing the electron and hole electric current densities J^ji to 
the corresponding electrochemical potential drops across 
the contacts of the device shown in Fig. [2] 



Je{L)^-^ViL)-V^ 
J,(L) = l|y(i)-Kf' 



tie{L)]} , (2.28a) 



e 



l^hiL)] 



Je(0) ^YvY'^ " e [^-^^^^ " 
MO) ^^{v^^ V{0) + - [m,.(0) + 



2r 



(2.28b) 
(2.28c) 

(2.28d) 



Here, t denotes the contact surface resistivity, i.e. Rc = 
x/W, where Rc is the electrical contact resistance and 
W is the sample width transverse to the temperature 
gradient. In Eqs. (|2.28ap - (j2.28dp . fiw denotes the single 
chemical potential level characterizing the particle-hole 
symmetric metal wire. Given the assumption of elec- 
troneutrality, remains fixed at its equilibrium value 
even in the presence of AT 7^ 0. 
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(see Fig. ig 



eVbi = -(Hw - Me) 





[ -Ih'wI 



eVbi = 



FIG. 3: This figure depicts contact ("built-in") potentials 
Vbi that arise through equilibration of graphene and the wire 
interconnect (voltmeter) depicted in Fig. [2] Subfigures (a) 
and (b) depict the electrochemical equilibration of electron- 
doped graphene and the metallic interconnect, while (c) and 
(d) show the same for the hole-doped case. All chemical po- 
tentials are measured from the position of the Dirac point; 
in this figure, we have assumed that fi^n > 0, i.e. the work 
function of the metal is less than the electron affinity of the 
graphene. The notion of a "negative hole chemical potential" 
—fivi for the carriers in the wire is easily understood via sub- 
figure (c): before electrochemical equilibration, the chemical 
potential difference — {—/j-w) drives holes in the graphene 
to "float up" toward the hole Fermi sea in the metal. Given 
the assumption of local electroneutrality throughout both the 
graphene and metal, the tunneling holes must recombine with 
electrons on the surface of the metal, leading to the accumu- 
lation of a dipole charge layer at the boundary. As a result, 
a static built-in voltage Vbi — —^{fJ-m + fJ-h) develops between 
the wire and graphene which precisely offsets the intrinsic 
chemical potential difference. Of course, this built-in voltage 
is not directly measurable with a voltmeter. 



In equilibrium, Ti ~ T2 ~ T; everywhere within the 
graphene slab, we have the conditions 



^e(a;) ^ -~lJ.h{x) ^1, V{x) 

Je{x) = Jh{x) = 0. 



0, 



Generically, the electron chemical potentials in the 
graphene slab (/^e = fi) and metal wire (/iu,) will differ; 
equivalently, ^.h — i—IJ^w) = —fi + l-iw 7^ 0. The interpreta- 
tion of fiui and — /iu) as the "electron" and "hole" chemical 
potentials in the metal is elaborated in Fig. [31 Assuming 
strong bulk screening in both materials, a dipole layer 
develops at the contacts, leading to the contact potential 



e 



In the nonequilibrium case, we therefore write 

= -{fl ~ fi^) + SV^, V° = -{fl - + SV^- 
e e 

(2.29) 

We now take the ideal limit of an infinitely resistive 
wire interconnect (Fig. [5]), and thus we require that the 
electric current vanish everywhere: 

J = 0. 



Using Eqs. ^JU^, (pUj) (j^l^ . and fTI^ . the bound- 
ary conditions in Eq. (|2.28p may then be recast as 

^uiL) =e\JniL), -Ai/(0) = e2rJ„(0), (2.30a) 
SV^ = V^{L), 6V° = 14(0), (2.30b) 

where we have introduced the electrochemical potential 
fluctuation 



1 



V<S>(x) EE V(x) + -^x). 

e 

Finally, the temperature fluctuation 'i>{x) satisfles 

AT 

vI'(O) = -^{L) = — . 
and $ were introduced in Eq. (|2.27p .] 



(2.31) 



(2.32) 



E. Solution to the linearized hydrodynamic 
equations 

Employing standard thermodynamic identities and lin- 
earizing in AT, we rewrite the hydrodynamic Eqs. (|2.6p - 
([2:81) . (|2T5l) . and ([2T6l) as the following set of 5 first 
order differential equations, valid to the lowest order in 



F- 

1 d^{x) 
T dx 



6VDa 



[N,Jnix)~2uN^ 



Vple 



dV,s,{x) /i d'${x) 
dx 



dx 



diJLi{x) 
dx 

dJn {x) 

dx 



+ 



2Do 



J„(x) -t- N,iU 



2D„ 



2n\( 



(Sf + E?) 



■/i/(x), 



Jn{.x)+N,U 



dx 



(2.33) 

(2.34) 

(2.35) 

(2.36) 
(2.37) 
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The various new parameters appearing in Eqs. (|2.33|) ~ 
([07)) are defined by 

= (Tee — CFeh, = ahh — CTeh, 

= {nlcThh + nlaee - 2nenh(7eh) ■ (2.38) 



We note that the kinetic coefficient sum 

EJ + S? = a^in, (2.39) 

the minimum conductivity at the Dirac point [Eqs. (|2.2ip 
and ((2:241) ]. 

When supplemented with Eqs. (|2.30ap and ((O^ . 
Eqs. (j2.33p - (|2.37p are easily solved. The results are 



4'(x) 



ErAT 



21q 



sinh 



f L-2x \ 

V 2/Q ) 



D„ L 



cosh 



(at? 



^min {-iV? (L-2X 



(2.40) 



ei 



j 2e 



e V L 



sinh (^^^' 



21q 



2Da L ^Qgj^ j 



(2.41) 



M/(a;) = ~ 1 ^ 



ST'TV.Si 2/q sinh(i-f) 



L 



cosh (2^) 

(2.42) 



in detail the inhomogeneity of the temperature and num- 
ber current density profiles implied by Eqs. p.40p and 
p.43p . General results for both n and the thermopower 
a are obtained and discussed in the final Sec. IIVI 



J„(x) 



f ) e^L 



cosh 



L-2x 
2lo 



cosh 



(2^3) 



AT 



3V ^ 



(2.43) 
(2.44) 



where Iq denotes the imbalance relaxation length, de- 
fined via 



1 _ nAQCTmm e2 



/2 



2D„ 



(2.45) 



In Eq. (CTD . 

Kao is the thermal conductivity in the limit 
of infinite imbalance relaxation, as given by Eq. (|2.22p . 
(See also Sec. IIIII below.) Finally, the parameters 91l 
and Si appearing in Eqs. p.40p ~ (|2.44p are defined as 



^L^r+ tanh — 

2ne^\olo \ 2ln 



(2.46) 



r^^tanh^ 



L \ , ^KLtT^in (37^)2 



^00 Tc^ 

" (2.47) 

In the next section, we specialize to the undoped case, 
and compute the thermal conductivity k in the limit of in- 
finite contact surface resistivity (r — > 00). We will discuss 



III. THERMAL CONDUCTIVITY AT THE 
DIRAC POINT 

We turn now to the calculation of the kinetic coef- 
ficients characterizing the thermoelectric transport. In 
this section, we concentrate upon the simplest case, that 
of zero doping. In equilibrium, the Dirac point is char- 
acterized by the conditions 



0, 



(3.1) 



As a result of the particle-hole symmetry implied by 
Eq. (|3.ip . the thermopower vanishes, a = 0. [This re- 
sult is demonstrated explicitly via Eq. (|4.2p in Sec. llVj . 

The thermal conductivity k obtains from the heat cur- 
rent [Eq. (|2.18p ]. Within linear response, 



J„ 



AT 

K 

L 

3-Pu-f 



O {ATY 



(3.2) 



In this section, we assume an ideal measurement of k, in 
which the contact electrical resistivity becomes arbitrar- 
ily large. 



00 



(3.3) 



[Sec Eqs. (|2.28aP " (j2.28dp .] Then, we use Eq. ((2:44)) and 
impose in addition the conditions hsted in Eq. (|3.ip upon 
all equilibrium thermodynamic variables, arriving at the 
result 



1 + L i-j^ ( L 

e^vpn'^l^l 21q V ^ 



where we have introduced 



cr hh 



2a. 



eh 



2 (a. 



(3.4) 



(3.5) 



As the so-defined n depends upon the sample length L, 
it is sometimes more natural to introduce the thermal 
conductance Gth, 



_ W 



(3.6) 



where W is the width of the graphene slab perpendicular 
to the applied thermal gradient. 

The thermal conductance in Eqs. (|3.6p and (|3.4p con- 
stitutes a primary result of this paper. The imbalance re- 
laxation due to non-electron and hole number-conserving 
inelastic collisions enters through the ratio L/Iq, where 
the length Iq was defined by Eq. (|2.45p . At the Dirac 
point, the latter simplifies to 



1 

13^ 



2n\( 



tjj' 



(3.7) 



The effects of quenched disorder are encoded in Eq. 
through the elastic mean free path /ei ■ 

Let us interpret our results. Consider first the clean 
limit, with /ei — oo. In this case, the physics is deter- 
mined by the ratio of the system size L to the length scale 
Iq [Eq. (j3.7p ]. For a short device, L ^ Iq, the thermal 
conductivity given by Eq. (|3.4p asymptotes to 



■■oo 




FIG. 4: The top graph (a log-log plot) depicts the qualitative 
form of the thermal "conductivity" k = LGtu/W [Eq. (|3.4|l ] 
versus sample length L for an undoped graphene strip pos- 
sessing both slow imbalance relaxation and weak quenched 
disorder. Three manners of functional behavior for k are de- 
marcated with dotted, solid, and dashed line segments; for 
each of these, the bottom plots show representative spatial 
profiles of the temperature T{x) - T = *(a;) [Eq. ([OSH and 
number current density Jn{x) [Eq. (|2.43|l ]. For system sizes 
Iq ^ L < loo (solid curves), k grows linearly with L, while 
the temperature profile is inhomogeneous; this regime is also 
characterized by the maximal particle flux J„ , as measured at 
the device center x — L/2. By contrast, in the short (L <C Iq) 
and long {L ^ loo) sample size limits, respectively designated 
by dotted and dashed curves in the figure, k saturates to Ko 
and Koo, respectively; here, the temperature profile asymp- 
totes to a linear gradient. The assumption of infinite contact 
resistivity [Eq. p.3p ] ensures that Jn{x) vanishes at a; = 0, L. 



o 



3^(3) 



L 



ITT- 



(3.8) 



The prefactor in this last equation obtains from the equi- 
librium pressure and density in Eq. p.4p . evaluated for 
the ideal relativistic quantum gas at zero doping, taking 
into account valley and spin degeneracies in graphene. 
In this short device limit, the temperature profile '^[x] 
[Eq. (|2.40p ] falls approximately linearly across the entire 
sample, while the number current density Jn is small ev- 
erywhere along the strip. 

In the opposite limit of a long device, L ^ Iq, the 
thermal conductance defined via Eq. (|3.6p approaches the 
L-independent constant 



r 

2<Q 



o 



cxp 



L 



(3.9) 



The temperature profile T{x) = T + '^{x) now consists of 
three regions: within a distance Iq of the sample bound- 
aries, the temperature drops approximately linearly; in 
between these boundary regions, T ^ T. A large car- 
rier number current J„(a;) [Eq. (|2.43p ] flows through the 
bulk of the sample, but pinches off to zero at a; = 
and X = L where, for AT > 0, generation and recom- 
bination processes, respectively, relax the accumulating 
population imbalance. As indicated by Eqs. (|2.43p and 
(|2.46p , if we relax the condition stipulating ideal contacts 
[Eq. (|3.3p ] by assuming a non-zero contact conductance 
density 1/r > 0, then J„(x) adopts a non-zero value at 
the sample edges; in this situation, electrons and holes 
that penetrate (escape) the graphene are generated (re- 
combined) in the contacts. 

For a long device {L » Iq) possessing, in addition, 
weak quenched disorder, the graphene sample behaves 
as three thermal resistors in series. The boundary resis- 
tance is dominated by the imbalance relaxation processes 
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as in Eq. (|3.9p . but there is now a finite temperature 
drop through the bulk. Disorder also introduces a second 
length scale loo ^ Iq into the denominator of Eq. (|3.4p . 
defined as 



(3.10) 



The scale loo emerges from Eq. p.4p when all thermo- 
dynamic variables in that equation are evaluated for the 
ideal quantum relativistic gas with zero charge density. 
For a graphene slab with L ^ loo, the bulk thermal resis- 
tance dominates, and the thermal conductivity k asymp- 
totes to the infinite imbalance relaxation value given by 
Eq. (|2.22p . which simplifies at the Dirac point to 



(3.11) 



For an appropriate definition of lei , Eq- (|3.1ip is the same 
result obtained in previous workfi^ which effectively as- 
sumed infinite imbalance relaxation. In Fig.Hl we sketch 
our results for the L-dependence of the temperature and 
number current profiles T{x) and Jn{x), as well as the 
thermal conductivity k = LGtu/W, for a sample satisfy- 
ing Ioo^^Iq- 

Note that in the clean limit lf.\ — > c», the result in 
Eq. (|3.1ip appears to suggest that the thermal conduc- 
tivity at zero doping diverges for any L. For any finite 
imbalance relaxation /g > 0, we have seen that this con- 
clusion is incorrect; instead, the response is inhomoge- 
neous (Fig. [4]), yielding a finite Gth for all L [Eqs. (|3.6p 
and ()3.4|1 ]. We observe that k in Eq. (j3.4|l depends only 
upon the combination or [Eq. ()3.5p ] of "intrinsic" kinetic 
coefficients; this is different from the combination that 
enters into the minimum conductivity at the Dirac point 
tTmin [Eqs. ([2?2T|) and ((2?24l) ]. Finally, it is important to 
stress that the limit given by Eq. (|3.1ip is appropriate 
only to the hydrodynamic transport regime, Tgi > T\n] in 
the opposite case, k is constrained by the "generahzed" 
Wiedemann-Franz lawp^ and a different result can be 
obtained. 



IV. THERMOELECTRIC TRANSPORT 
COEFFICIENTS: ARBITRARY DOPING 

In this section, we compute the thermoelectric re- 
sponse of graphene at non-zero doping. Combining 
Eqs. (|3.2p and (|2.44p . the general expression for the ther- 
mal conductivity is 



(3P)^„ 

fe2 



(4.1) 



The thermoelectric power a was identified in Sec. |TT] 
as the ratio expressed in Eq. (|2.26p . where AX^j ~ 
5V^ — 5V^ is the voltage drop across the "wire" intercon- 
nect serving as our voltmeter in the experiment sketched 



in Fig. [2l The boundary conditions in Eq. (|2.30bp demon- 
strate that a is determined by the electrochemical po- 
tential drop AVi ~ V<s>{L) — Vb(0) across the graphene. 
Eq. (|2.4ip then implies that 

A14. 



AT 



(4.2) 



^,VS,NAX^-^2_l, 

Te 2D„ L \2l 



The parameters S^, Nj, Tfj'^ , and D„ appearing in 
Eqs. (|4ll) and gj) were defined in Eqs. ([235]) . jMH), 
and ((Mil), above. 

The general expressions in Eqs. (|4.ip and (j4.2p consti- 
tute the primary results of this paper. We now specialize 
these results to the short (L ^ Iq) and long {L ^ Iq) 
device limits. 



1. Short device L <^Iq 

For a device shorter than the imbalance relaxation 
length Iq [Eq. ^T^], Eqs. (jil^ and (ji?^ asymptote 
to 



L 



(3-P)V^in 



Olf 



3P_ 



(4.3) 



AveNij + Lp 



Te 



Iq 
L 



(3P)2 



o 



(4.4) 



These results indicate that the limit L Iq gives not one 
but several regimes for the thermoelectric response in the 
general case; the demarcation lines between these depend 
upon L and the contact surface resistivity r, as well as 
the doping and the extent of disorder in the sample. 

We analyze Eqs. (|4.3p and (|4.4p . neglecting the influ- 
ence of disorder for simplicity (^ei — *■ oo) . We focus upon 
the non-dcgcncratc case of high temperatures and low 
doping, ksT 3> A maximum of three behavioral 

regimes are possible for both k and a, and these are 
accessed sequentially with increasing L. We define two 
length scales mediated by the contact resistivity r, 



l' 



II 



(4.5) 



where n and p arc the total number and charge densities, 
respectively [Eq. (|2.12p ]. In Eq. (|4.5p . fTmin denotes the 
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minimum dc conductivity at the Dirac point in the hy- 
drodynamic regime, as defined by Eqs. (|2.2ip and (|2.24p . 
Clearly we have Z^"^' ^ 4*^' in the non-degenerate regime. 
The interpretation of If''^'' is as follows. At the scale lf\ 
the electrical conductance of the graphene sample at zero 
doping is of order the contact conductance, since 



W<T„ 



1 



(4.6) 



where Rc = r/VF, and W is the sample width trans- 
verse to the current flow. By contrast, is the scale 
at which the thermal conductance Gth for clean, non- 
degenerate graphene in the infinite imbalance relaxation 
limit becomes of order the electrical contact resistance: 



Gth.oo {L — /J. 



klf 



1 



(4.7) 



where k,oo is given by Eqs. (|2.22[) and (|2.21[) in the limit 
^ei — > oo. 

For a device with L <^ lf\ we can set L = ev- 
erywhere in both Eqs. (|4.3p and (|4.4|) . The resulting 
expressions for k and a are exactly those obtained for 
graphene possessing zero imbalance relaxation, Aq = 
in Eqs. (|2.16p and (|2.45p . as measured through "ideal" 
(electrically insulating, thermally conducting) contacts. 
These zero imbalance relaxation (shortest device) expres- 
sions for K and a are completely different from the infi- 
nite imbalance relaxation (longest device) results and 
ctoo given by Eqs. (|2.22p and (|2.23p . above. For exam- 
ple, using the definitions of N,i and N"^ from Eq. (|2.38p . 
it is clear from Eq. (|4.4p that the thermopower in this 
regime vanishes smoothly as the sample charge density is 
tuned through the Dirac point, even for a perfectly clean 
sample. By contrast, Eq. (|2.23p gives aoo ~ s/pva the 
limit Zei ^ oo, which exhibits a simple pole at p — Q (s 
is the entropy density). Measurements of the thermo- 
electric response in both the short and long device lim- 
its for graphene in the hydrodynamic regime should al- 
low extraction of the independent diagonal (uee) and off- 
diagonal {a eh) tensor coefficients defined by Eq. (|2.16[) . 

In a longer device with intermediate L, If^ <C L ^ /t*^' , 
the terms proportional to L in the numerators of both 
Eqs. (|4.3p and (|4.4p dominate the response. Similar to 
the intermediate regime of the thermal conductivity of 
undoped, disordered graphene discussed in Sec. IIIII and 
illustrated in Fig. 31 both k, and a rise linearly with L. 
For much longer devices with L ^ 4''' , the term propor- 
tional to L in the denominator of both Eqs. (|4.3p and 
(|4.4p dominates, so that k and a plateau to their infi- 
nite imbalance relaxation limits {Koc>,aoo), transcribed 
above in Eqs. (|2.22p and (|2.23p . In other words, when 
the contact electrical conductance becomes larger than 
the infinite imbalance relaxation limit of the graphene 
thermal conductance Gth.oo (times e'^/kgT) [Eq. (|4.7p ]. 
we find that k — s- Kqo. 

Observation of the described crossovers requires that 
the imbalance relaxation length Iq exceeds one or both 



of if'-' <^ if'-' . In the non-degenerate regime of the carrier 
plasma, we obtain an order of magnitude estimate for Iq 
by approximating Eq. (|2.45p as 



(4.8) 



[See also Eq. (|3.7p . Sec. Illlj . For a fixed sample size L, 
we can determine the temperatures T = T'r'^''^' below 
which the crossovers at L if-' and L if^ become 
observable within the short device (L <C Iq) regime. As 
it is expected to vary only weakly with decreasing tem- 
perature, the smaller scale If^ may be approximated as 
a constant. Eq. (|4.8p then implies that the crossover at 
L if'-' occurs for temperatures 



2^ < = 



hvp 1 



(4.9) 



By contrast, the f-dependence of ^t'"' is partially deter- 
mined by the conditions of the experiment. If for example 
the charge density p = e(ne — nn) is held constant as the 
temperature is varied, then Eq. (|4.8p suggests that the 
crossover &i L ^ If^ should be observable in the short 
device regime only for temperatures 



T < ri' 




(4.10) 



These equations hold only for relatively clean, non- 
degenerate graphene; the average chemical potential and 
temperature must satisfy \p,\/kB <C T. In addition, if 
we take Aq ^ 1 and dmin 4e^//i (see the discussion in 
the next subsection, below), then Eq. (|4.10p holds only 
for Zei ^ Iqi^n/ p)'^ ■ By contrast, for /ei ^ Iq, only 
one crossover within the L <^ Iq regime is possible, at 



2. Numbers for the short device (L <^ Iq) regime 

Experimentally, vp ~ 10^ m/s, while a^in ~ Ae^/h:^^ 
We have not calculated the imbalance relaxation parame- 
ter Aq in Eqs. (|4.8p - (|4.10p . which requires a quantum ki- 
netic equation treatment that incorporates three-particle 
collisions and/or impurity- assisted recombination. We 
have argued that, due to the absence of two-particle 
mechanisms (see Fig. [T] and the concomitant discussion 
in Scc.|T|, imbalance relaxation should be a slow process 
in the hydrodynamic regime. In what follows, wc take 
the conservative estimate Aq ~ 1. 

We consider first the crossover at L ^ If '' . For a con- 
tact resistance oi Rc = 100 and a relative carrier den- 
sity ofp/e = 10i°/cm2, Eq. gH) gives T,'"' ^ 400K. The 
average (relative) chemical potential is p./kB ^ 2QK for 
the assumed density^ while If^ - Zq(TT') ~ 0.02 ^m. 
Longer crossover lengths if^ and lower temperatures 
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Tj'^' can be obtained for larger contact resistances Rc, 
but lower carrier densities are required to preserve the 
condition of non-degeneracy for the carrier plasma. A 
contact resistance of Rc = lOfcfi and a relative car- 
rier density of p/e ~ 10^/cm^ gives T^'"^' ^ AK, with 
If' - IgiT^"'') - 2 i^im. In this case, p/kB 0.2K. 

We now turn to the crossover at L ~ li'^K For 
a contact resistance Rc = 1^1 and a relative density 
p/e = 10"/cm2, Eq. ^iJU^i gives T^*"' - 700A'. The 
chemical potential is /i/fc^ ~ lOOii' for the assumed 
density^ while 4"' ~ ^g(rr') 0.01 ^ini. In order 
to observe deviations from the infinite imbalance relax- 
ation limit [Eqs. (|2.2ip - (|2.23p ] in the short device regime, 
the sample length L < so these conditions would 
seem to require an impractically short device. Longer 
devices meeting the required constraints arc possible for 
more resistive contacts and lower carrier densities. For a 
contact resistance of Rc = 100 fl and a relative density 
p/e = lOVcm^, one finds T,'"' ~ 7K, p/kg IK, and 
it'' ~ Iq (Tr ) ~ 1 pm. Contact resistance can be con- 
trolled in principle through the incorporation of a highly 
insulating spacer layer of varying thickness between the 
graphene and the contact metal. 



3. Long device Iq 

In the opposite limit where the sample length L ex- 
ceeds the imbalance relaxation length Iq, Eqs. (|4.ip and 
(|4.2p simplify as follows: 



8rK^ Te^Nj 
(37')V„i„[4D„r+/QCr„,i„] 



O [exp {-L/Iq)] , 



(4.11) 



Kp 



ST'tTmin Te 

'Iq\ en 4riV,(S? - E^') 



0[exp(-L/ZQ)]. (4.12) 



The denominator of Eq. (|4.1ip introduces yet another 
length scale into the problem. 



Iq + V. 



(4.13) 



where if''^' were introduced in Eq. (|4.5p . On the second 
line of Eq. (|4.13p . we have used Eq. (|2.22p . neglected 
the effects of disorder for simplicity, approximated ^ 
~ (Tmin and D„ ^ c^ip, and we have assigned to all 
thermodynamic variables their values at the Dirac point. 
In the limit L > ?oo, Eqs. (|nT|) and ((iTT^ show that 
the thermal transport coefficients asymptote toward their 



infinite imbalance relaxation limits, k — > and a{K,) — > 

The criterion for the existence of a length-dependent 
crossover in the behavior of k and a within the L ^ Iq 
regime is as follows. For Iq :p if , Eq. (|4.13p gives 



(4.14) 



The scale is also the crossover scale to the same (ef- 
fective) infinite imbalance relaxation limit, as obtained 
in the opposite regime L <^ Iq. Thus, the location of the 
crossover to infinite imbalance relaxation behavior rela- 
tive to L ^ Iq depends upon the unspecified ratio iQ/li'"', 
consistent with the previous discussion. 

In the opposite limit Iq <C llf ' , a crossover in the L ^ 
Iq regime definitely occurs at 



L 



Iq\- 



(4.15) 



with n [p) the total number (charge) density. In this 
case, an intermediate regime exists for Iq <^ L <^loo, in 
which K and a grow linearly with L. Only for L ^ loo do 
the infinite imbalance relaxation limits for these kinetic 
coefficients emerge. 

We stress once again that Eqs. (j4.3p . (|4.4p . (|4.1ip . and 
(j4.12p hold only for the case of hydrodynamic. interpar- 
ticle collision- mediated transport. In the opposite case of 
"disorder- limited" transport, where the elastic scattering 
rate exceeds the inelastic rate due to interparticle colli- 
sions (rei ^ Tin), K and a are slaved to the electrical con- 
ductivity through the "generalized" Wiedemann-Franz 
law and Mott relation, respectively.— Further discus- 
sion on the distinction between interaction and disorder- 
limited transport in graphene can be found in Ref. [ill . 

Finally, we comment upon the physics of the thermo- 
electric transport within the L ^ loo regime. From 
Eq. (|2.42p . we note that, in the long device limit, 
the imbalance chemical potential pi{x) [introduced in 
Eq. ()2.13p ] is exponentially suppressed between boundary 
layers of size Iq. The electron- hole population imbalance 
is therefore confined to the boundary regions. By con- 
trast, Eq. (|2.43p shows that the number current J„ that 
fiows through the bulk of the sample decays only linearly 
with increasing L, for fixed AT: 



J„(L/2) 



2Nj 



AT 



L knT' 



(4.16) 



where we have used Eq. (|4.ip and approximated '--^ Sj*. 
For L ^ Iq, ti cc L aX the Dirac point for clean graphene 
[Eq. (|3.4p with /el oo]; at zero doping J„ therefore 
saturates to a finite, non-zero value as the system size 
diverges, consistent with the picture of the central region 
as a "perfectly conducting thermal wire." By contrast, 
the 1 /L decay of Eq. (|4.16p away from the Dirac point is 
consistent with the finite thermal drop implied by in 
Eqs. (|iTT|l and ^TF^. 
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V. CONCLUSION 

111 summary, we have demonstrated that thermoelec- 
tric transport in graphene within the hydrodynamic 
regime exhibits a range of behaviors when the finite 
rate of carrier imbalance relaxation is taken into ac- 
count. Since the relativistic spectrum of clean graphene 
is non-decaying, the lowest order two-particle recombina- 
tion and generation processes are kinematically forbid- 
den, suggesting that the imbalance relaxation lifetime 
TQ might significantly exceed other intrinsic graphene 
timescales. 

The essential transport physics in the hydrodynamic 
regime is encoded by four intrinsic parameters: these 
are the minimum conductivity at the Dirac point tTmin, 
the off-diagonal (or "drag" ) conductivity a eh , the imbal- 
ance relaxation length Iq , and the elastic mean free path 
le\. Of these, (Tmin, cFeh, and Iq are mediated entirely by 
intercarrier collisions. The parameters cTmin and le\ can 
be obtained from measurement of the bulk conductiv- 
ity (at variable doping), or the combined measurement 
of the electronic thermal conductivity k(L) and the ther- 
mopower a{L), in the limit of a long device with L ^ Iq. 
The drag conductivity aeh can be extracted from a mea- 
surement of either k or a in the opposite, short device 
limit L Iq. For a sample with L > Iq, a local probe 
of either the electronic temperature or electrochemical 
potential profiles should allow determination of Iq , since 
these are predicted to be inhomogeneous, with boundary 
layers of size Iq confined near the device terminals. 

We have given general formulae for both k and a at 
arbitrary doping and device size L, incorporating the ef- 
fects of non-ideal contacts. Non- ideal contacts allow ex- 
change of carriers with the graphene, providing an alter- 
nate route for imbalance relaxation. We have explicated 
the various crossover regimes that separate the zero im- 
balance relaxation (short device) and infinite relaxation 
(long device) limiting behaviors. 

Ill this paper, we have neglected the effects of phonons 
in graphene. Both acoustic and optical phonons can in- 
fiuence the electronic thermal conductivity contribution 
K and the thermopower a, through inelastic elcctron- 
phonon scattering. Specifically, real electron-phonon col- 
lisions may modify (i) the imbalance relaxation rate due 
to electron-hole pair to phonon conversion processes, (ii) 
the inhomogeneous electronic temperature profile, due 
to energy exchange with the phonon bath, and (iii) the 
thermoelectric power through phonon drag. By contrast, 
virtual electron-phonon interactions are strongly irrele- 
vant, and the concomitant renormalization effects may 
be typically neglected. 



For temperatures less than Swph/fcs TOOK, all opti- 
cal modes are frozen-out Untethered ("free floating") 
graphene supports linearly-dispersing acoustic phonons 
within the transverse (TA) and longitudinal (LA) in- 
plane modes, as well as quadratically-dispersing phonons 
in an out-of-plane (ZA) mode. Under tension imposed 
by external contacts or surface adhesion to a substrate, 
however, the ZA dispersion also becomes linearii 

Because the acoustic phonon velocities^ iiph ~ 10^ m/s 
^ vp, the electron-hole pair creation and annihilation 
processes 

+ ^ ph 

are kinematically forbidden. Therefore the electron- 
phonon scattering does not contribute to the imbalance 
relaxation, at least to lowest order. For graphene in the 
non-degenerate regime with /cgT ^ ftwph, the in-plane 
acoustic phonon contribution to the inelastic electron 
lifetime (due to electron and hole number-conserving pro- 
cesses) Te_ph oc {kBT)~'^, by dimensional analysis. The as- 
sociated electron-phonon relaxation length ^e-ph = VpTe-ph 
can be estimated with the Boltzmann transport resulli^ 

, 4(fa>F)3p„,«g^ 

^-'^ ^ {ksTDf ' ^^^^ 

where pm 7.6 x 10"^ kg/m^ denotes the 2D mass den- 
sity of graphene, Wph ~ 2 x 10^ m/s is the phonon veloc- 
ity for the LA mode, and D ~ 19 eV is the deformation 
potential.— Using these parameters, Zg-ph ^ 80 /im at 
T = lOOK and 0.8 cm at T = lOK. For devices with 
sample dimensions L,W < Ze-ph, phonons may be ne- 
glected. In larger devices, the loss of carrier energy to 
the phonon bath on scales longer than l^-pu becomes im- 
portant, so that the parameter /g-ph will enter e.g. into 
the temperature profile across the device. In addition, 
phonon drag effects can become important in sufficiently 
clean samples with L ^ Ze-ph- 
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